van der Waals interactions between graphitic nanowiggles 
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The van der Waals interactions between two parallel graphitic nanowiggles (GNWs) are calcu- 
lated using the coupled dipole method (CDM). The CDM is an efficient and accurate approach 
to determine such interactions explicitly by taking into account the discrete atomic structure and 
many-body effect. Our findings show that the van der Waals forces vary from attraction to repulsion 
as nanoribbons move along their lengths with respect to each other. This feature leads to a number 
of stable and unstable positions of the system during the movement process. These positions can be 
tuned by changing the length of GNW. Moreover, the influence of the thermal effect on the van der 
Waals interactions is also extensively investigated. This work would give good direction for both 
future theoretical and experimental studies. 

PACS numbers: 



I. INTRODUCTION 



The van der Waals interaction, arising from the quan- 
tum fluctuations within bodies, has received enthusias- 
tic attention in recent years due to its important role 
from condensed matter physics [TH3] to biological sys- 
tems |4, 5J. At nanoscale, this force is much more domi- 
nant than other forces such as the electromagnetic force 
[5J. In the field of material science, the van der Waals 
force mainly causes the stiction and friction problems 
which directly impact behaviors and physical properties 
of nanodevices. The understanding of such interactions 
has provided the cornerstone towards next-generation de- 
vices and advancement fundamental science. 

There are several main ways to estimate the van 
der Waals interactions. The well-known method imple- 
mented in molecular dynamics (MD) simulation codes 
is based on the Lennard-Jones (LJ) model (7J-[5]. This 
calculation, however, considers only the pairwise interac- 
tion between two atoms and ignores the effects induced 
by other atoms. Another approach is built up from the 
Lifshitz theory [TU1 - [T2"] . This theory allows us to calcu- 
late the van der Waals and Casimir force between two 
parallel semi-infinite plates using dielectric functions of 
objects as a function of the imaginary frequency. The 
interaction between two objects with arbitrary geome- 
tries such as sphere-sphere and cylinder-cylinder systems 
can be calculated by the proximity force approximation 
(PFA) [THl E3] . A different heuristic approach for study- 
ing dispersion interactions derives from the multipole ex- 
pansion [14] . In the framework of this method, the van 
der Waals energy between two spherical bodies with the 
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centre-to-centre distance R are expressed by 



where C 6 , C 8 , and C w ,... are the van der Waals co- 
efficients that can be obtained by using the dynamical 
polarizability of each object. To apply Eq.Q effectively, 
the separation of two spheres is required to be larger 
than a certain minimum value. In practical complex sys- 
tems such protein, DNA and graphene nanoribbons, the 
anisotropy and surface roughness must be taken into ac- 
count in the van der Waals interactions. As a result, 
the long-range interaction at short distances estimated 
by these two above methods may not accurate. In this 
paper, we use the CDM to calculate the van der Waals 
interaction energy. This method is based on the many- 
body effects and has been used broadly in a wide range 
of research areas [T5HT7] . The CDM is capable of dealing 
with systems containing a large number of atoms and the 
effects caused by the structure orientations. 

It is well-known that graphene is an extraordinary ma- 
terial due to its fascinating properties and possible ap- 
plications 18J. Graphene also provides state-of-the-art 
classes of material such as carbon nanotubes (CNTs), 
bulky balls, graphene nanoribbons (GNRs) and GNWs 
(see Fig|T]). Recently, both experimentalists [TU^O] and 
theoreticians [2~TH2"5] have been attracted by GNWs. The 
finite-size effects and the modification of the edge in the 
chevron nanostructurc influence significantly its prop- 
erties and lead to be better performance than pristine 
graphene. The increase of electron-hole confinements 
strengthens the exciton binding energy over that of regu- 
lar GNRs [23]. One can take advantage of this feature by 
designing optical sensors to be used in biosystems |24j . 
Moreover, it was shown that the chevron-like GNRs have 
large values of the figure of merit and can be useful in 
thermoelectric devices [22] . In this work, we focus on the 
van der Waals interaction in chevron-like GNRs. 
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L = 9 
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For short separations kRij <C 1, the retardation effects 
can be ignored. In this range, the highest order terms of 
1 / R dominate over others. The magnetic contributions in 
Eq.|2]), therefore, is neglectable. Eq.([2]) can be rewritten 



as pi 
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in which x^, yij and are components along x, y and 
z axis of Rij, respectively. It is important to note that 
Ty = as i = j. ai is the atomic polarizability of dipole 
i modelled by the Drude model [22] 



(6) 



where aoi is the static polarizability and ujQi is the char- 
acteristic frequency of the ith dipole. In our work, the 
atomic polarizability is assumed to be isotropic. Now it 
is easy to obtain the 3N coupled equations [2"T] 



FIG. 1: (Color online) The atomic arrangements of a graphene 
nanowiggle. The grey circles denote carbon atoms while the 
white circles refer to hydrogen atoms. The length of GNW is 
specified by the number of L. 



II. THEORETICAL BACKGROUND 

This computational theory starts by considering a sys- 
tem as set of N discrete polarizable atoms interacting to 
each other via electromagnetic fluctuations. A polariza- 
tion of the dipole at position r j is induced by [25] 
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Solving Eq.Q provides 3N values of u>. The van der 
Waals interactions of the entire system at if is ex- 
pressed by 
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where fj,Q is the magnetic permeability of the vacuum, 
nij is the magnetic dipole moment of the jth dipole, the 
free-space Green's function is given 



The retarded dispersion energy between object 1 and 
2 are calculated by 
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where k — lo/c is is the wave number, c is the speed of 
light, I is the 3x3 identity matrix, = r\; — Tj is the 
distance vector between dipole i and j, R^ = ILy/.Rjj, 
and the vector is expressed 



g« = 
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where Uf w , Uf w and U vdw are the van der Waals 
energy of only object 1, object 2 and the whole sys- 
tem including both of them, respectively. In some cases, 
the tensor given by Eq.([5| leads to a "polariza- 
tion catastrophe" problem. To avoid this issue, it is 
necessary to introduce the distance-dependent function 
a (aoiaoj) 1 ^ 6 and modify the coupled dipole inter- 
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where a — 1.662 is the screening factor. If R^ > Sy then 



Vij = 1, else if Rij 



< Sij then Vij 



Rij I Si 



3 



III. NUMERICAL RESULTS AND DISCUSSION 

In this section, we investigate the van der Waals in- 
teractions between two chevron-like GNRs in vacuum. 
The used parameters for atomic polarizabilities j37J EH] 
are u c = 1.85 x 10 16 rad/s, w# = 1.41 x 10 16 rad/s, 
a c = 0.85 A 3 , and a H = 0.25 A 3 . Each sheet is re- 
laxed by molecular dynamics (MD) simulation ReaxFF 
and then we combine these two ribbons to create a static 
system. 
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FIG. 2: (Color online) The van der Waals interactions of two 
GNWs versus the movement of top nanoribbon along the y 
axis at K. The dashed line illustrates the fitting function 
of the van der Waals energy around y=0 at z = 4 A. 

Figure [2^ shows the dispersion energy between two 
parallel identical chevron-like GNRs with L = 9 at sev- 
eral certain separation distances z along the z axis using 
the CDM method. As the separation z increases at the 
same y, this interaction decreases significantly. Drifting 
the top nanoribbon along the y axis causes an oscillation 
of the retarded dispersion energy. This oscillatory-like 
feature reveals that the van der Waals force changes from 
attraction to repulsion and vice versa. At z = 4 A, the 
oscillation is so strong, and we can see the stable equilib- 
rium positions of the system at local minimum energies 
corresponding to y = 0, 18 and 34 A. Obviously, the 
most stable equilibrium state is at y = for all distance 



z, other positions have feeble barriers and can be called 
the metastable positions. At long separation distances, 
the oscillatory-like properties of energy reduces and al- 
most disappears at z = 8 A. Because the influence of 
atomic structures, particularly stacking, is significantly 
weakened when the separation distance z is beyond the 
C-C bond length. 

The results in Fig|2^ also present that the oscillatory 
trend between two nearest peaks is approximately 18 
A. This distance corresponds to the salient geometri- 
cal period of the chevron-like GNRs depicted in Figjl] 
The interpretation for this peculiarity is that the x — y 
plane projection of the top nanoribbon match with that 
of the bottom one. The oscillatory feature can ob- 
tained using the Lennard- Jones (LJ) 12-6 potential = 
-A/Rfj+B/Rjf shown in Figj^D. Here A and B are the 
Hamaker constants. The parameters Acc — 15.2 eVA 6 
and B C c = 24 x 10 3 eVA 12 representing the C-C inter- 
action in the graphene-graphene system was introduced 
in Ref.[31 . The parameters for the C-H and H-H inter- 
actions A'ch = 11-4 eVA 6 , B CH = 11.77 x 10 3 eVA 12 , 
A HH = 8 eVA 6 and B HH = 5.05 x 10 3 eVA 12 were cal- 
culated from the MD simulation such as LAMMPS code. 
Although the L-J model is the simplest method to es- 
timate the interaction between two neutral atoms and 
has been widely used in numerous MD codes, the pair- 
wise approach does not provide highly accurate results 
in comparison with the CDM due to the underestima- 
tion of the many-body effects. Our findings and previous 
study [32 show the difference of the magnitude of the 
van der Waals interactions derived from the L-J model 
and the CDM method. One, however, can observe the 
same phenomenon by these two calculations. 

Using fitting function, the approximate analytical ex- 
pression of the van der Waals energy Ufa around z = 4 
A and y = is given in Fig. [2^,. The form of this fitting 
function indicates that we can consider the small move- 
ment of top GNR along the y axis as harmonic motion. 
The frequency of this oscillation can be easily calculated 
by 
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where m is the mass of the top nanoribbon. One can 
estimate m = 3y/3na 2 p g /2, in which a = 1.42 A is the 
bond length between two carbon atoms, p g = 7.6 x 10~ 7 
kg/m 2 is the surface mass density of graphene [33], and 
n = 132 is the number of hexagon in GNWs. Substi- 
tuting all parameters to the expression of m, we obtain 
m — 5.255 x 10~ 24 kg and the frequency vq ~ 2.42 x 10 11 
Hz. Another method to calculate m is the sum of mass 
of individual atoms in the ribbon |34j and this approach 
leads to m = 7.76 x 10~ 24 kg. However, this estima- 
tion does not take into account the reduction of total 
mass due to the formation of the binding energy when 
bringing isolated atoms closer together. At retarded dis- 
tances, the van der Waals interactions dominates other 
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forces. The motion of the top flake is controlled only by 
this long range interactions. 

In the next step, we study the effect of GNW length on 
the van der Waals energy. The length of bottom GNW 
remains unchanged = 9 while the length of top GNW 
varys with L t = 6, 7, 8 and 9. As can be seen in Fig. 
[3J the dispersion interactions between two GNWs at the 
same separation z = 4 A the reduction of L t leads to 
the unbalance of top nanoribbons at y = 0. The stable 
equilibrium positions are approximately 8 A and 4 A at 
L t — 7 and 6, respectively. For L t = 8, the local minima 
of the energy interactions at y = and 12.5 A have the 
same energy value, however, the energy barriers around 
y is smaller than those of y = 12.5 A. Therefore, y = 
12.5 A would be the most stable equilibrium position of 
this system. Our findings indicate that the finite size 
structure and the edge state have an immense impact 
on the van der Waals energy. Furthermore, there are 
a number of repulsive-attractive transitions of the van 
der Waals force at y < 36 A. At larger distances, only 
attractive van der Waals force exists. 
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FIG. 4: (Color online) The van der Waals interactions of two 
GNWs versus the movement of top nanoribbon along the y 
axis at zero and room temperature when the center-center 
distance of two objects in z direction is 4 A. 
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FIG. 3: (Color online) The van der Waals interactions of two 
GNWs versus the movement of top nanoribbon along the y 
axis. 

It is remarkable that one can modify Eq.([8| to esti- 
mate the van der Waals energy of the whole at certain 
temperature T 

where ks is the Boltzmann constant. At room tem- 
perature, ksT ss 0.025 eV and all eigenvalues provide 
huji/k B T » 1. As a result, Eq.([8]) can be applicable to 
the calculation at non-zero temperature in this case, so 
17$*% = U^ dw in terms of formula. The interesting prob- 
lem is that the thermal effect influences the roughness 
of the nanowiggle. The different rough surfaces result in 
different energy interactions. 



To understand the impact of the thermal effect on the 
van der Waals energy, we calculate the interactions sim- 
ilar way to previous calculations, however, each GNR is 
relaxed separately by ReaxFF program at 300 K. It can 
be seen in Fig. |4j the van der Waals energy interactions in 
the system including the small GNW are more sensitive 
than systems possessing the longer chevron-like GNR. It 
is well-known that the bond in edge is weaker compared 
to the rest area on the flake. Interestingly, temperature 
induces the bend of the GNW edge. This deformation 
generates the energy difference in the system under differ- 
ent thermal conditions. For short nanoribbons, the edge 
effect plays an important role in the interaction between 
two objects. This effect, however, has less contribution 
as the length of GNW is increased. For GNW L t = 9, 
the dispersion energy is nearly unchanged at y < 10 A 
as temperatures vary. While one can observe the large 
discrepancy of such interaction as T = K and T = 300 
K. At large y, the energy difference at different temper- 
ature in both cases L t = 6 and L t = 9 tend to disappear. 
This result implies that the edge A is far apart enough 
to ignore the influence of edge atomic coordinates on the 
van der Waals interactions. 



IV. CONCLUSIONS 

We have discussed the van der Waals interaction be- 
tween two chevron-like GNRs by means of CDM. The in- 
fluence of stacking and edge effect causes the oscillatory- 
like feature of this energy as the top flake moves parallel 
to its length direction. Local minima and maxima of en- 
ergy curve is shifted when the length of a nanoribbon is 
varied. This coincides with the change of the most sta- 
ble position of the system. Heating up systems induces 
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the variation of edge atomic positions. At retarded dis- 
tances, this structure change provides the difference of 
van der Waals energy at different temperatures. As the 
center-center separation of two GNWs increases, the en- 
ergy discrepancy becomes smaller. 
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